Weighted Assortative And Disassortative Networks Model 
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Real-world networks process structured connections since they have non-trivial vertex degree 
correlation and clustering. Here we propose a toy model of structure formation in real-world weighted 
network. In our model, a network evolves by topological growth as well as by weight change. In 
addition, we introduce the weighted assortativity coefficient, which generalizes the assortativity 
coefficient of a topological network, to measure the tendency of having a high- weighted link between 
two vertices of similar degrees. Network generated by our model exhibits scale-free behavior with a 
tunable exponent. Besides, a few non-trivial features found in real-world networks are reproduced 
by varying the parameter ruling the speed of weight evolution. Most importantly, by studying the 
weighted assortativity coefficient, we found that both topologically assortative and disassortative 
networks generated by our model are in fact weighted assortative. 

PACS numbers: 89.65.Gh, 05.65.+b, 05.70.Fh, 89.75.-k 
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I. INTRODUCTION 

There is a wave of interest to study complex networks 
by statistical physical means EL!.2, 3]. The area of stud- 
ies include the Internet [1, the World-Wide Web 
Q, scientific collaboration networks (SCN) biolog- 
ical networks [l(| E| and the world-wide airport network 
(WAN) [T2|. In all these cases, one can naturally identify 
the subject under study with a graph. For instance, an 
airport in a WAN can be represented by a vertex; and 
there is a link between two vertices if and only if there 
is a direct flight between the corresponding airports. Al- 
though these complex networks are drawn from vastly 
different systems in the real world, they exhibit the fol- 
lowing universal properties: 

1. Scale-free behavior of degree distribution: Let P(k) 
be the probability that any vertex is connected to 
k other vertices (i.e. the vertex is of degree k). In 
many real-world networks, P(k) ~ fc -7 with 2 
7^30. 
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2. Clustering property: Clustering of a vertex Ci gives 
the probability that two nearest neighbors are con- 
nected with each other. And the average clustering 
coefficient C = J2i c i/N , where N is the number of 
vertices in the network, measures the global den- 
sity of interconnected vertex triplets in the network. 
Real-world networks in general exhibit higher clus- 
tering coefficient than those in random networks 

3. Small-world property: A network is said to have 
small-world property if the average shortest path 



length between two vertices scales at most loga- 
rithmically with N for fixed mean degree [3, EH ■ 

4. Correlations of degree: An assortative (A disassor- 
tative) network tends to connect vertices with sim- 
ilar (dissimilar) degrees. Previous studies found 
that social networks tend to be assortative while 
technological and biological networks are generally 
disassortative [lq . 

Various network models have been proposed to sim- 
ulate or explain the properties found in complex net- 
works in the real world. For instance, the Erdos-Renyi 
random graph model [13, EH generates static networks 
that exhibit small-world property. Nonetheless, the de- 
gree distribution of the network generated is Poissonian 
rather than a power-law. Later on, Barabasi and Albert 
(BA) extended the Erdos-Renyi random graph model by 
evolving the network using a linear preferential attach- 
ment mechanism (l3j. The ideas behind the BA model 
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are that most real-world networks tend to growth with 
time, and a new vertex tends to connect to pre-existing 
high-degree vertices. Sen further incorporated non-linear 
growth property of the number of links to the BA model 
by demanding the number of links of the vertex added 
at time t to be t e for some positive constant 9 • Both 
the BA model and the Sen model generate networks with 
scale-free topology; but they failed to reproduce the de- 
gree correlations and clustering properties of real-world 
networks. 

Structural organization of a network can be character- 
ized by degree-dependent average nearest-neighbors de- 
gree k nn (k) [J], the degree-dependent average clustering 
coefficient C{k) [2(| and assortativity coefficient r [2lj |. 
Pan et al. proposed the generalized local-world models 
[22| in which a subset of vertices is randomly chosen every 
time step. The newly added vertex can only make con- 
nection to vertices in this subset. Moreover, Liu et al. 
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proposed the self-learning mutual selection model [23| in 
which every vertex has some probability to self evolve in 
each time step. These two models successfully reproduce 
the k nn (k), C{k) and r observed in real- world networks. 
However, these quantities only depend on the topological 
structure of a network. The relative importance between 
different links is not taken into consideration. To pro- 
vide a better description of a network, the intensity of 
interaction among vertices should be taken into account. 
One may define the weight of a link u)y as the intensity 
of interaction between vertices i and j, and the vertex 
strength Sj as the sum of weight of the links connected 
to i: 



jer(i) 



(1) 



II. WEIGHT EVOLUTION MODEL 

The model starts with a small number Nq of fully 
connected undirected vertices. The weight of all the 
No(Nq — l)/2 initial links are set to wo. In each time 
step, the network evolves under the following rules: 

1. Topological growth: A new vertex is added to the 
network. We assume that the number of links of 
the new vertex is an increasing function of the net- 
work size. Specifically, a vertex with [pt e ] links are 
added to the network at time t for some fixed con- 
stants p > and < 9 < 1. (Here [x] denotes the 
value of x rounded to the nearest integer.) These 
[pt e ] links are all of weight u>o and are randomly 
connected to the existing vertices according to the 
strength preferential probability Aj, which is de- 
fined as & m 
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where T(i) denotes the nearest neighbors of i. For in- 
stance, the weight in WAN refers to the number of 
seats available on the direct flight connections between 
the airports i and j. And s,; represents the total traf- 
fic going through the airport i. To probe the weighted 
networks' architecture, a set of weighted quantities are 
introduced, such as the weighted degree-dependent aver- 
age clustering coefficient C w (k) and the weighted degree- 
dependent average nearest-neighbor degree k™ n (k) [2J|. 
However, no one has proposed a quantity that can di- 
rectly measure the correlations of degree with inclusion 
of weight so far. Here we introduce the weighted assor- 
tativity coefficient r w to measure the tendency of having 
a high-weighted link between two vertices of similar de- 
grees. In addition, we propose a model that incorporates 
the essential features of strength preferential attachment, 
nonlinear growth of number of links, and weight evolu- 
tion of existing links. We call this the Weight Evolution 
Model. The rule of weight evolution are based on the no- 
tion that "the rich always gets richer". In other words, 
high-weighted link has higher probability to evolve. Be- 
sides, inspired by the work of Dorogovtsev and Mendes 
[25l | . existing links in the graph can be removed. By al- 
tering the speed of weight evolution, our model can gen- 
erate both assortative and disassortative networks with 
scale-free behavior and non-trivial clustering properties. 
Most importantly, both topologically assortative and dis- 
assortative networks generated by our model are in fact 
weighted assortative. 

In Section [H] we define the Weight Evolution Model. 
We report our analytical calculations and numerical re- 
sults on the probability distributions of strength of vertex 
and weight of link in Section IIII1 Then in Section IIVI we 
define the weighted assortativity coefficient and compare 
the organizational structure found in real-world networks 
with our numerical simulations. Finally, we give a brief 
summary in Section fVl 
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2. Evolution of weight: The link between vertices i 
and j is chosen to evolve with probability Qij pro- 
portional to its weight, namely 



(3) 



J2a<b W ab ' 

The weight will update according to the equation 

Wij — ► Wij + sgn(g) , (4) 

where 



1 if q > 
sgn(q) = ^ —1 if q < 
if q = 



(5) 



At time step t, this process will repeat 
times and a link is removed if the weight equals 
to 0. (In other words, the sign of q determines 
the addition or subtraction of weight of a link; and 
the magnitude of q specifics the rate of the addi- 
tion/subtraction per time step.) 

Note that p, q and 9 are the tunable parameters in this 
model. 

Clearly, our model generalizes the BA model [l3| and 
the Sen model In particular, our model is reduced 
to the BA model by setting 9 = Q, q = 0, Wo = 1, and 
reduced to the Sen model by setting < 6 < 1, p = 1, 
q = and Wq = 1. 

We now justify the validity of our rules. Note that in 
most real-world networks, new vertices are introduced 
one at a time; and they tend to connect to the pre- 
existing high strength vertices. In the example of WAN, 
new airport tends to establish flights to existing airports 
with heavy traffic. Our topological growth process is 
used to capture this kind of network growth. Moreover, 
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our weight evolution mechanism models the trend that 
the link with higher weight have a higher probability to 
evolve. For instance, in SCN [24j], a vertex represent a 
scientist and a link is present if two scientists have co- 
authored at least one paper. The weight between two 
scientists is higher if they have higher number of co- 
authored papers. Naturally, if the weight between two 
scientists is high, the probability for them to collaborate 
again is also high. 



III. STRENGTH AND WEIGHT 

In our model, the processes of topological growth and 
weight evolution are Markovian. Thus, the evolution of 
our network can be calculated accurately by mean field 
approximation. Appendix [X] reports the mean field cal- 
culation of the probability distribution of vertex strength 
P(s) and probability distribution of weight of link P(w). 



We found that P(s) 



ls = l- 
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and P(w) 



2(w p + q) 
{0 - l)w p - 2q 



w p 



' , where 



1)9 



(6) 



(J) 



To check the validity of our analytical solution, we per- 
formed numerical simulations on our network model for 
different values of q and 8. Since p, q and wo govern the 
relative speed between new link addition and weight evo- 
lution, one can fix two of these three parameters without 
losing any generality. In this paper, we fix p — 1 and 
wo = 10 unless otherwise stated. Besides, we set the 
number of initial vertices Nq to 5; and we have checked 
that the networks generated using different values of p 
and No exhibit similar behaviors as long as Nq <C N. Us- 
ing this set of parameters, we found that statistical data 
collected from the network equilibrate after t > 5 000. 
Hence, we collect all our statistical data at t = 10 000 
and all the data reported here have been averaged over 
50 independent runs. Our model successfully reproduces 
the scale-free behavior of the probability distributions of 
strength and weight with a tunable exponent that de- 
pends on the microscopic mechanism ruling the weight 
evolution. We found in Appendix [B] that the numerical 
simulation results agree with the mean field approxima- 
tion. 



IV. CLUSTERING AND CORRELATIONS 

In this section, we first discuss the quantities used to 
characterize the clustering and correlations behaviors of 



a network. Then we introduce the weighted assortativity 
coefficient r w . Finally, we report the numerical simula- 
tion results of our model. 

A. Quantities Characterizing A Weighted Network 

We overload the notation by labeling the vertex added 
to the network in time step i also by i. The clustering of 
i is defined as 



3,h 



ihO-jh 



(8) 



where ki is the degree of i and 



_ 1 if a link exist between i and j , , Q -, 
i if there is no link between i and j . * ' 



If vertex i has less than two neighbors, Ci is set to 0. 
Clustering Cj measures the local cohesiveness of vertex 
i while the average clustering coefficient C — ^ i Ci/N 
measures the global density of interconnected triples in 
the network. Organizational structure of network can be 
further studied via the degree-dependent average cluster- 
ing coefficient C(fc), 



C(k) = 



NP(k) 



E < 

i:ki—k 



(10) 



which is the mean clustering for vertices with degree k. 

To have a better understanding of the organizational 
structure of weighted networks, the weighted clustering 
cf is introduced and defined as 1241 



i{ki 



E 

3,h 



W lh ) 



~(lijCL%}\ Q>jh 



(ii) 



If vertex i tends to form triples with other vertices by 
its' high-weighted links, then cf > a meaning that the 
topological clustering of i underrates the cohesiveness of 
i. The weighted average clustering coefficient C w and 
weighted degree-dependent average clustering coefficient 
C w (k) are defined as the average of cf over all vertices 
and over all vertices with degree k, respectively [2~jj . 

The degree correlation is another important informa- 
tion of the network. Recall that the network is said to be 
assortative if the vertices tend to connect to other ver- 
tices which have similar (dissimilar) properties. Newman 
introduced the assortativity coefficient [16L [2l| 
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to measure the degree correlation between linked vertices, 
where F{<j>) denotes the set of the two vertices connected 
by the 0th link and M is the total number of links in 
the network. This measure r is positive (negative) for 
assortative (disassortative) networks; and r = for a 
random graph. Note that r is independent of the weight 
of each link of a network. 

To probe the degree correlation of the network, one 
may also study the average nearest-neighbors degree, 
which is defined as [H 

k n n,i — ~T~ ^ t kj- (13) 

i ier(i) 

The degree-dependent average nearest-neighbors degree 
k nn (k) is the mean of fc nn ,i restricted to the class of de- 
gree k vertices. In an assortative (disassortative) net- 
work, vertices with high degree tend to connect to other 
vertices with high (low) degree, thus k nn (k) would be 
an increasing (decreasing) function of k. In real-world 
weighted networks, high degree vertices could connect to 
small degree vertices with low weight, while connect to 

I 



where Wi is the weight of the 0th link, F{<j>) is the set of 
the two vertices connected by the 0th link and H is the 
total weight of all links in the network. Just like r, r w lies 
between —1 and 1. Moreover, r w is positive for weighted 
assortative networks, while negative for weighted disas- 
sortative networks. If r w > r, a high- weighted link tend 
to connect two similar degree vertices together. In fact, 
r w reduces to r if all weights in the network are equal. 
Furthermore, r w = r = for a random graph. Fig. [1] 
illustrates that the values of r and r w of a network can 
differ greatly. The network in Fig. Q] consists of two com- 
plete graphs of three vertices that are joined together by 
a highly weighted link. Since the five out of seven of 
the links in this network connect two vertices of different 
degree together, the assortativity coefficient r of this net- 
work is negative. (In fact r = —1/6.) On the other hand, 
the weighted assortativity coefficient r w of this network is 
positive. This can be understood as follows. The link be- 
tween the two degree three vertices in this network carries 



I 

high degree vertices with high weight. For instance, in 
WAN, the high degree airport i could have a lot of direct 
flight to another high degree airport j, while have less 
number of flight to a low degree airport h. In this case, 
k nn (k) will underestimate the tendency for having heavy 
traffic between two similar degree airports. To handle 
this problem, Barrat et al. proposed the weighted aver- 
age nearest-neighbor degree [24| 

If the weighted degree-dependent nearest-neighbors de- 
gree k™ n (k) is an increasing function of fc, similar degree 
vertices tend to link together. Besides, the weights of 
these links tend to be high. Consequently, the network 
is weighted assortative. 

Although one can figure out that the network is assor- 
tative or disassortative if k™ n (k) is increasing or decreas- 
ing with k, a quantity directly describing the weighted 
assortativity is needed. Thus we introduce the weighted 
assortativity coefficient 



(15) 

I 

the most weight. So after coarse-graining, the network is 
similar to the one making up of a single link connecting 
two degree one vertices together plus another four iso- 
lated vertices. And clearly, the coarse-grained network is 
assortative. Therefore, when the weight of link are taken 
into account, the network in Fig.Q]ought to be assortative 
rather than disassortative. (In fact r w — 2/3.) In this 
respect, the weighted assortativity coefficient r w better 
determines the assortativity of networks with weighted 
links than the assortativity coefficient r. 



B. Simulation Results 

As shown in Fig. [21 our model can generate assortative 
networks (r > 0), such as those found in social networks, 
when q 0. Our model can also generate disassorta- 
tive networks (r < 0), such as those found in techno- 
logical networks, when q 3> 0. Surprisingly, r w > 
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w=10 
w=l 



r= -0.167 r w = 0.333 



FIG. 1: A network whose r < but r w > 0. 
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FIG. 2: Assortativity coefficient r and weighted assortativity 
coefficient r w versus q for various 6. The statistics taken right 
at time t = 10 000. And we have averaged over 50 independent 
runs in all the data reported here. All error bars are < 10 -2 . 



for both assortative and disassortative networks provided 
that |g| > 7. This observation indicates that the topolog- 
ical assortativity coefficient underrates the contributions 
of high- weighted links between similar degree vertices. 
Our finding shows that studying the assortativity of a 
network by considering network topology alone [H, 
does not give a complete picture of the organizational 
structure of the network. 

The difference between r w and r (Fig. [3]) can be under- 
stood qualitatively by considering the evolution mecha- 
nisms of our model. For ij > 0, a newly added vertex is 
the one with low degree and it tends to connect to ex- 
isting high-strength and high-degree vertices. This pro- 
cess leads to topologically disassortative behavior. Mean- 
while, a link introduced at an early time, which connects 
two old vertices having similar degrees together, tends to 
be high- weighted as q > 0. This leads to the observed 
weighted assortative behavior. On the other hand, as- 
sortative networks emerge at q <C 0. It is because the 
strength of old vertices are low. Hence, there is a high 
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FIG. 3: The comparison between r and r™ for 8 = 0.3. 
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FIG. 4: fe nn (fc) and kZx(k) for various q for 9 = 0.3. 



chance for two young vertices having similar degrees to 
connect. As 9 increases, the strength of old vertices de- 
crease while the new vertices enter the network with a 
higher strength. Thus, for <; C 0, the value of \q\ required 
for the emergence of assortative behavior decreases with 
9. This is precisely what we have observed in Fig. O 

As shown in Fig. [U k nn (k) exhibits decreasing power- 
law behavior for g>0, indicating that the networks are 
disassortative. For g«0, k nn (k) is a gently increasing 
function of k. It is worth noting that in some empiri- 
cal studies of real- world networks [24j . both k nn (k) and 
k™ n (k) show assortative behavior with fc™ n (fc) > k nD (k). 
As shown in Fig. our model successfully reproduces 
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FIG. 5: Comparison of k nn (k) and fc™ n (fc) for (a) q — 10 and 
(b) q = -9 for 6 = 0.3. 
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creasing power law for q ^> 0. Moreover, C w (k) is larger 
than C(k), especially at large fc, which is a feature found 
in real- world networks [HI ■ This result indicates that the 
topological clustering coefficient C(k) underestimate the 
cohesiveness of the network. It is because C(k) cannot 
tell us whether large degree vertices tend to form inter- 
connected triples with its' high-weighted links. 



V. CONCLUSIONS 

In summary, we have introduced the Weight Evolution 
Model which couples dynamical evolution of weight with 
topological network growth. We have also introduced 
the weighted assortativity coefficient r w to characterize 
the degree correlation of a weighted network. Our model 
reproduced many features found in real-world networks, 
such as scale- free behavior, assortative and disassortative 
behavior, and two types of clustering structures. Most 
importantly, both topologically disassortative and assor- 
tative networks are regarded as weighted assortative if 
the weight of links is taken into account, indicating that 
the topological assortativity coefficient r underrates the 
contributions of high-weighted links between similar de- 
gree vertices. This result may give us some idea for why 
social networks are found to be assortative, while techno- 
logical networks and biological networks show an oppo- 
site behavior in previous studies [la. |2l| . It is instructive 
to investigate if all real-world networks can be regarded 
as weighted assortative. 
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FIG. 6: C(k) and C w (k) for various q for 6 = 0.3. 



these important features. Remarkably, r w < r while 
k™ n (k) > k nn (k) for q <c 0. This points out that it is 
difficult to compare the weighted assortativity and topo- 
logical assortativity by studying the difference between 
k™ n (k) and k nu (k). 

Two types of C{k) are found in real- world networks 
[2fl [28j . In the first type, C(k) does not exhibit strong 
dependency on k. While in the second type, C(k) shows 
a decreasing power law spectrum. Fig. [S] clearly shows 
that our model can reproduce both kinds of features. In 
particular, C(fc) is flat for g < and it becomes a de- 



APPENDIX A: ANALYTICAL CALCULATIONS 
OF STRENGTH AND WEIGHT 

Recall that initially the network is a complete graph 
of iVo vertices. And in each time step a new vertex is 
added. There are two ways to alter the strength of an 
existing vertex i: (a) i is connected to the newly added 
vertex j with probability Aj given by Eq. or (b) the 
weight of the link between i and one of its neighboring 
vertices evolves according to Eqs. ([3])-([5]). 

By mean field approximation and by treating all dis- 
crete variables as continuous, the strength Si of the vertex 
added at time i satisfies 



dsj(t) 
dt 



= w pt 



i(t) 



<lt' 



er(») w n 



(t) 



J2a<b W ab(t) 
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w pt 



««(*) 



2qt b 



Sjjt) 

J2i s i(t) 



(Al) 



for t > i. Note that the first and second terms repre- 
sent the strength change due to topological growth and 
evolution of weight, respectively. 

Since the total strength of the network increases 
2(wop + q)t B approximately in each time step, 



i 



2{w p + q)t e dt 



2{w p + q) 4e+1 



1 



(A2) 



Putting Eq. ![A"2j) into Eq. ([AT]) and using the initial con- 
dition Sj(f — i) — Wo[pi e ], we conclude that 



Si{t) 



w pi 



w pt 



■ \ (6+l)(w p+2q)/2(w p+q) 



,\ [2g-(6l-l)«)op]/2(ti)of>+«) 



(A3) 



Since one vertex is added to the network in each time 
step, the network size Nfzt. The probability distribution 
of vertex strength can be computed by 



i(t))dj , 



(A4) 



where S(x) is the Dirac delta function. Combined with 
Eq. (|A3|) . we conclude that in the infinite network size 
limit, the mean field approximation give P(s) ~ s~ ls for 
min(si(t)) < s < max(sj(t)), where 



7s 



1 



2QoP + g) 
1 - l)w p - 2q 



(A5) 



The behavior of P(s) can be classified according to the 
value of q: 



1. If q > q c = (9 - l)io p/2, Eq. jA3J) tells us that the 
degree of existing vertices are generally larger than 
the degree of the newly added vertex, i.e. Si(t) > 
s n = wopt e for all i. Thus, 



P(s) 



for s < s„ 
s _7s for s > s n 



(A6) 



In particular, the network evolves purely by topo- 
logical growth for q = 0. And in this case, we find 
that 7 S = (6 — 3)/(0— 1) which is independent of 
p. Moreover, 7 S — > 2 as q — > oo. 

2. If (7 = (7c, Sj(f) = s n which is the same for all i. In 
this case, P(s) is a delta function. In contrast, P(s) 
is a delta function only when 9 = 1 (i.e., when the 
network is fully connected) in the Sen model [l9[ . 



3. If q < q c , Si(t) < s n by Eq. ()A3j) . In other words, 
the newly added vertex have the highest strength. 
Using a similar argument as in the case of q > q c , 




FIG. 7: Probability distribution of strength P(s) for 6 = 0.3. 
A power-law function P(s) ~ s ' is obtained in the range of 
(a) s < s n if q < q c = —3.5 and (b) s > s n if q > q c - The 
dashed line corresponds to s = s„ = 148. The insert compares 
the values of 7 obtained from data fitting (filled circles) with 
the analytical solution 7 S given by Eq. (1A5[) (solid line). 



we found that the strength distribution P(s) fol- 
lows 



P(s) 



s 7s for s < s n , 
for s > s n . 



(A7) 



The evolution of weight can again be computed using 
mean field approximation. Therefore, we have 



dwjj (t) 
dt 



qt" 



Wij(t) 



a (9 + l)qwij(t) 

J2a<b W ab(t) " (W p + q)t 



(A8) 



for t > max(«,j). Clearly, Wij(t) satisfies the initial con- 
dition Wij(t = max(i, j)) — Wq. Hence, we find 



Wij(t) w wo 



-I (e+\)q/(w p+q) 



max(z,j) 



(A9) 



In addition, the weight distribution P(w), which gives 
the probability that a link with weight w, obeys P(w) ~ 
f or mhi(wij) < w < max(wij), where 



lw = 1 



w Q p 



(0 + l)q 



(A10) 
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FIG. 8: Probability distribution of weight P(w) ~ w with 
(a) q < and (b) q > 0. The insert compares the values of -y w 
obtained from data fitting (filled circles) and the mean field 
value given by Eq. (|A10[) (solid line) for 6 = 0.3. 



APPENDIX B: NUMERICAL RESULTS OF 
PROBABILITY DISTRIBUTION 



Fig. reports P(s) obtained from numerical simula- 
tions for different values of q. It also compares the value 
of exponent obtained from fitting the numerically simu- 
lated data with the value predicted by analytical calcu- 
lation. P(s) follows a power-law function for g < g c or 
q ^> q c with exponent agree well with the value given by 
Eq. (|A5|) . Note that P(s) is non-zero for s > s n when 
q < q c (Fig.[7Ja)) and for s < s n when q > q c (Fig. [7(b)). 
This discrepancy is due to the fact that a small propor- 
tion of vertices evolves differently from the prediction of 
mean-field approximation. However, this value dimin- 
ishes to zero rapidly and thus our analytical calculation 
remains valid. The weight probability distribution P(w) 
for various q is displayed in Fig. [8] along with the com- 
parison between the fitted values of the exponent with 
the analytical predictions. For q < 0, lUy < wo for all 
links, so the power-law function is obtained for w < wq 
(Fig. [8(a)). For q > 0, the power-law function is found 
to be in the range of w > w (Fig. [S^b)). 
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